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ABSTRACT 

We present the first study of the far- infrared (FIR) properties of high redshift, radio- 
selected Ultraluminous Infrared Galaxies (ULIRGs) using deep observations obtained 
with the Spectral and Photometric Imaging Receiver (SPIRE) from the Herschel 
Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 
850/im fluxes from submillimetre-luminous ~10 mJy sources (SCUBA galaxies) to 
~1.5 mJy from stacked SCUBA non-detections, thus likely representing a complete 
distribution of ULIRG spectral energy distributions (SEDs). From Keck spectroscopic 
surveys in the Lockman-North field we identified a sample of 31 submillimetre galaxies 
(SMGs) and 37 submillimetre- faint, optically-faint radio galaxies (OFRGs), all with 
radio-inferred IR luminosities > I0 12 L©. These galaxies were cross-identified with 
SPIRE 250, 350 and 500/im catalogs based on fluxes extracted at 24/im positions in 
the SWIRE survey, yielding a sample of more than half of the galaxies well detected 
in at least two of the SPIRE bandpasses. By fitting greybody dust models to the 
SPIRE photometry together with SCUBA 850/im measurements (for OFRGs, only 
850/im upper limits), we infer dust temperatures and far- infrared luminosities. The 
OFRGs detected by SPIRE have median (T dust ) = 41 ± 5 K and the SMGs have 
(Tdust) = 34 ± 5 K, both in reasonable agreement with previous (pre- Herschel) esti- 
mates, reaffirming that the local FIR/radio correlation holds (at least for this subset 
of high- z ULIRGs) at high redshift (we measure (<7ir) = 2.43 ±0.21 using Sir, derived 
from greybody fit coupled with a powerlaw extrapolation to the 24/iin). Our obser- 
vations firstly confirm that a substantial fraction of OFRGs exhibit large infrared 
luminosities corresponding to SFRs of ~ 400 M Q /yr. The SPIRE observations sec- 
ondly confirm the higher dust temperatures for these OFRGs than similarly selected 
SMGs, consistent with early predictions of the submm-faint radio populations. Our 
observations also clearly confirm the large infrared luminosities of most SMGs selected 
with Sg50/im > 5 mJy and radio and strong 24/im detections, corresponding to SFRs 
of - 700 M Q /yr. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: infrared — 
galaxies: starbursts 



1 INTRODUCTION 

Submillimetre Galaxies (SMGs - Smail, Ivison & Blain 
1997) contribute significantly to the rapid buildup of stellar 
mass in the Universe at z ~ 2. SMGs have a typical 
redshift of z ~ 2.2 (Chapman et al. 2005; Wardlow et al. 
2010), are massive systems (M, ~ lO lo ~ 11 Af , Swinbank 
et al. 2004; Greve et al. 2005; Tacconi et al. 2006) and 
are diverse in the extent and dynamics of their molecular 
gas reservoirs (e.g., Tacconi et al. 2008; Bothwell et al. 
2010; Ivison et al. 2010a,2010b). Interferometric observa- 
tions of SMGs' CO molecular gas suggest that the most 
luminous SMGs S (850/im) > 5 mJy are merging sytems 
(Engel et al. 2010) with high star formation efficiencies 
compared to typical galaxies of similar mass (Daddi et 
al. 2008). However, their selection at 850/im is inherently 
biased towards c o lder-dust ULIRG s, particularly at z>l 
l|Eales et al. 2000l : iBlain et al. 20041 ). In particular, since 
submm observations probe the blackbody emission of dust 
in the Rayleigh- Jeans regime, they are anti-correlated with 
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dust temperature (Ssso oc T^^JJ for a given infrared lu- 
minosity, and galaxies with warmer dust can fall below the 
dete ction limit of current submm instruments. Recent work 
(e.g. IChapman et al. 2004al : ICasev et al. 2009) has demon- 
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1 since L oc M d T 5 - 5 while the Rayleigh-Jeans flux Srj oc M d T. 



strated that 850/im- faint, high-redshift ULIRGs exist and 
may contribute significantly to the cosmic star formation 
rate density at its peak. These Optically Faint-Radio Galax- 
ies (OFRGs) are defined as radio sources having inferred 
ULIRG luminosities, with starburst (SB) or hybrid SB-AGN 
spectral features in the UV, and having 2.5-<r limits on their 
submm fluxes which are consistent with them being fainter 
than 5 mJy at 850 /im. They have a comoving volume den- 
sity (i.e., ~ 10~ 5 Mpc~ 3 at 1 < z < 3, Chapman et al. 2001, 
2004), stellar masses and radio sizes comparable to SMGs, 
and some have a dust temperature of ~52 K (Casey et 
al. 2009, 2010a, 2010b). Studies of other i nfrared-luminous 
galaxy populations bo th pre-Herschel (see | E>ev et al. 20081 : 



lYounger et al. 2009| ; Bussmann et al. 20091) and post 



Herschel (see lOliver et al. 20ld: iMagdis et al. 20T0I : 



iRoseboom et al. 20ld : IChanial et al. 20101 ') present even 
more evidence for diverse populations of luminous, dusty 
starbursts at z ^ 1. 

However, sparse infrared data, particularly in the 50- 
500/im wavelength range have limited the interpretation of 
the SMGs and OFRGs. Many of their fundamental proper- 
ties still rely on indirect measurements. Direct determina- 
tions of SMG and OFRG dust temperatures are limited and 
have only been done using either a single rest-farlR point 
(SHARC-II/CSO observed 350/rni - Kovacs et al. 2006,2010; 
Coppin et al. 2008) or with the most luminous examples 
of the population using the highly confused BLAST beam 
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(250/xm, 350/im & 500/im Ivison et al. 2010c, Dunlop et al. 
2010, Chapin et al. 2010, Casey et al. 2010b). 

Observations of the rest-frame far-infrared emission of 
high-redshift galaxies with Herschel (Pilbratt et al. 2010) 
provide for the first time direct measurements of typical 
SMGs and OFRGs total infrared luminosities and dust tem- 
perature^, which are still highly debated (e.g., Swinbank 
et al. 2008; Dave et al. 2008). Importantly, they allow for 
tests of earlier measurements based on extrapolation from 
wavelengths either shortward or longerward of the farlR 
peak. Theoretical simulations of galaxy evolution continue 
to have difficulties in accounting for the inferred luminosi- 
ties/star formation rates and number counts (Baugh et al. 
2005; Dave et al. 2009; Narayanan et al. 2009). Open ques- 
tions remain as to whether these luminosities have been 
overestimated or whether for instance the IMF is signi- 
cantly more top heavy than in the local Universe (Baugh 
et al. 2005). In this paper, we assess the Herschel-SPIRE 
properties of these two ULIRG populations selected using 
both submillimetre (submm) and radio wavelengths. Calcu- 
lations assume a flat, ACDM cosmology with J1a = 0.7 and 
H = 71kms" 1 Mpc- 1 . 



2 SAMPLE & OBSERVATIONS 

The HerME^ Science Demonstration Phase (SDP) obser- 
vations are detailed in Table 1 of Oliver et al. (2010). We use 
deep Spectral and Photometric Imaging Receiver (SPIRE) 
(Griffin et al. 2010) 250, 350, 500/rni observations of SMGs 
and OFRGs lying in the Lockman-East (lOh 51m, +57) re- 
gion (a subset of the larger Lockman-SWIRE field observa- 
tions) obtained as part of this guaranteed time key program 
in a key survey field having exquisite multi-wavelength ancil- 
liary data. The radio/submm sources lie within a sub-field 
covered by a deep VLA 1.4GHz radio pointing, within a 
region of 15' x 15' (Ibar et al. 2009). Flux densities are mea- 
sured using a PSF fitting technique based on cross-identifed 
(XID) prior sources positions detected at 24/im (Roseboom 
et al. 2010), which accounts for blending of SPIRE source 
fluxes by simultaneous fitting to multiple sources. All sources 
in our fields are essentially point sources at the resolution 
of Herschel. The association between the Spitzer-IRAC & 
MIPS and SPIRE sources is thus facilitated by the use of 
these priors. The quality of the SPIRE 250, 350, 500Atm 
catalogs has been estimated using Monte Carlo simulations 
(photometry error, completeness and contamination as a 
function of the flux density). In the Lockman-East field 
our observations have average la XID source extraction er- 
rors of 3.6 mjy, 4.2mJy and 4.8mJy at 250, 350, 500/im 
respectively, and the confusion noise is quoted as 5.8, 6.3 
and 6.8mJy/beam rms respectively by Nguyen et al. (2010), 
which we apply in quadrature to the flux errors. 

Some SMGs and OFRGs have previously been stud- 
ied spectroscopically in Lockman-east in several works (re- 
spectively, Chapman et al. 2003a, 2005; Ivison et al. 2005 
and Chapman et al. 2004; Casey et al. 2009,2010). The 

2 Although for the hotter and lower redshift examples in our 
sample, PACS observations on the Wien side of the spectrum are 
required to truly constrain the dust temperature 

3 hermes.sussex.ac.uk, Oliver et al. 2010 (in prep.) 




Observed Wavelength (/im) 

Figure 1. Average SEDs of SMGs (top) and OFRGs (bottom), 
after redshifting each template fit to the rest-frame, and interpo- 
lating observed data points at the average redshift of the sample. 
Circles present our SPIRE measurements, SCUBA 850/^m and 
VLA 1.4GHz data. Error bars show the dispersion in the sample, 
which is comparable to the individual flux errors combined with 
confusion noise. Red, green, and blue lines (segment, dashed, dot- 
ted) present the Dale & Helou SED templates from 25-45K. An 
average SMG template from Pope et al. (2008) is also overlaid as 
solid black line. 



parent sample for this study consists of a sample of 376 
radio sources with spectroscopic followup in the Lockman- 
east field. For this initial foray into the SPIRE properties 
of luminous, high-z radio sources, we impose a > 5a 24/im 
detection threshold (from the SWIRE catalogue - Lonsdale 
et al. 2004) in order to extract the SPIRE fluxes at that 
position. This removes some of the confusion issues arising 
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in the full radio source sample, although it brings some bias 
to the study. We adopt a catalogue based on a 24^im prior, 
as opposed to a radio prior, as we have uniform 24/xm cov- 
erage across these areas and a clear understanding of the 
effect of this input list on the data. Future contributions 
will included a broader radio prior catalog to further explore 
these populations. This catalog consists of 140 sources. We 
further restrict this sample to those with radio luminosi- 



ergs/s/Hz, corresponding to an equivalent 
1 1 N assuming the radio- 



ties > 4.5 x W 

ULIRG far-IR luminosity (> lO^L© 



farlR correlation (e.g. Condon et al. 1991), resulting in a 
sample of 68 sources. 

The final sample consists of 31 securely identied SMGs 
and 37 OFRGs with IR luminosities > 10 12 Lq, and which 
are cross-matched with our SPIRE 250, 350, 500/im multi- 
wavelength catalogs based on sources extraction using prior 
detections at 24^m. 19 of the SMGs and 21 of the OFRGs 
in our sample are detected in at least two of the SPIRE 
bandpasses, and form the sample we focus on in this work. 
The SMGs have an average (deboosted using the Coppin 
et al. 2006 counts and method) S850nm=6mJy while the 
OFRGs have average S850M™ = l-5mJy (which we have not 
attempted to deboost here), together likely representing a 
complete distribution of ULIRG SEDs. In figure 1 we plot 
the average SEDs of SMGs and OFRGs which is constructed 
by redshifting each template fit to the rest-frame, and inter- 
polating observed data points at the average redshift of the 
sample. The dispersion in the average SED is comparable 
to the individual flux errors combined with confusion noise. 
Figure 1 also compares the Dale & Helou SED templates 
ranging from 25-45K, and an average SMG template from 
Pope et al. (2008) . Figure 3 shows the redshift distribution. 
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Figure 2. The dust temperature-luminosity relation. Filled 
blue/empty red circles are from our SMGs/OFRGs Lockman-East 
sample. Temperature uncertainties range from 3-7 K (depending 
on how well the peak is sampled) while luminosity uncertainties 
are ~0.2 dex (given the power law extrapolation down the Wien- 
side model fit). The cyan squares represent results for SMGs ex- 
trapolated by Chapman et al. (2005) from radio and submm data. 
The Chapman et al. (2003b) derivation of the median and in- 
terquartile range of the Td-LiR relation observed locally is shown 
by solid and dashed-dotted lines, linearly extrapolated to 10 13 Lq . 



3 RESULTS 

We fit the SPIRE (250^m, 350^m, and 500/im), and SCUBA 
(850^im) (only the upper limit in the case of the OFRGs) 
flux densities to FIR dust models (see Figures 1 & 2). The 
model assumes a modified blackbody emission curve, in the 
optically thin approximation, with a single dust tempera- 
ture: 



exp(hv / kTdust) - 1 

where S v , the flux density, is a function of rest frequency 
v, the emissivity /?, dust temperature Td, and FIR luminos- 
ity Lfir (which governs the normalization of the function). 
Their total infrared luminosities (Lir [8 — lOOO^im]) are in- 
ferred from these best fit using the far-infrared luminosity 
denition (Lfir[40 — 120/im]) given by Helou et al. (1988) and 
a color-correction term (Dale et al. 2001, Lir = 1.91 x Lfir). 
The model fixes emissivity to ft = 1.5. The model has Td 
and Lir as free parameters. While the luminosities here are 
clearly an approximation to the true Lir, we leave a more 
detailed estimation to the availability of shorter wavelength 
Herschel data. We choose to make this model rigid since the 
flux errors are known to be correlated which reduces the 
effective number of degrees of freedom. 

This single dust temperature characterisation provides 
a good description of the SMGs and OFRGs far-infrared 
SED, however for 7 SMGs and 9 OFRGs, this single dust 



temperature model yields high \ 2 values. All these galax- 
ies appear either to be the less luminous examples, exhibit 
far-infrared colours possibly suggestive of hotter systems, or 
may simply lie in confused regions. However, given the ex- 
treme dust masses (~ 10 9 Mq) of SMGs, one should expect 
optical depths to matter, especially in the shorter wave- 
length SPIRE bands. Assuming typical emission scales of 
~2 kpc, we expect optical depths of 1-2 in the 250/xm and 
350/im bands for a z ~ 2 SMG (see Kovacs et al. 2010). 
While the effect is still moderate such that the optically thin 
approximation is crudely correct, the high \ 2 found for some 
of the galaxies may be an indication of the shortcoming of 
the single-T, optically thin model. 

Our results appear to be relatively insensitive to ft and 
to the single dust component model; with ft = 2, the differ- 
ences in Td are ATd ± 3 K. In addition, the multiple dust 
components model of Dale & Helou (2002) yields rest-frame 
infrared colors (S60Mm/Sioc>Mm)i or equivalently Td , which 
are in reasonable agreement with those inferred from our 
greybody analysis. In figure 1, we also overplot a compos- 
ite SMG spectrum, from Pope et al. (2008), normalized to 
24^im flux density. 

The Td-LiR plane forms an important diagnostic for 
assessing the properties of luminous IR galaxies and their 
spectral shapes. In Figure 2, the locations of our SMGs and 
OFRGs on the Td-LiR plane are shown. These estimates 
are amongst the first direct observational measurements 
of the dust temperatures and the infrared luminosities of 
OFRGs. Previous studies with SHARC-II and BLAST have 
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Figure 3. Submm flux densities as function of the redshift and 
Far-Infrared luminosity, revealing both the redshift distribution of 
our sample and the range in 850/im fluxes probed for two samples 
which have similar average radio fluxes. Symbols are same as in 
Fig. 2 except the OFRGs for which we have only upper limits 
(flux+2cr). Solid and dashed lines show the linear fit to the Ssso/jm 
- Lir relation for SMGs and the lcr envelope. 



attempted similar analyses mainly on SMGs with caveats 
presented in the Introduction (Kovacs et al. 2006, 2010, 
Coppin et al. 2008; Ivison et al. 2010c, Dunlop et al. 2010, 
Chapin et al. 2010, Casey et al. 2010, Amblard et al. 2010). 
Our observations suggest that high redshift ULIRGs show a 
wide range of dust temperatures. At intermediate infrared 
luminosities of ULIRGs (Lir ~ 3 — 5 x 10 12 Lq), the dust 
temperature dispersion observed in our sample might sug- 
gest a higher Td-LiR scatter than deduced from local IRAS 
galaxies by Chapman et al. (2003) and Chapin et al. (2009). 
OFRGs tend to be biased towards hotter dust tempera- 
tures, with median (Td) = 41 ± 5 K (quoting standard de- 
viations of the sample) and (Lir) = 3.8 x 10 12 Lq, while 
SMGs appear to have lower dust temperatures with median 
(T d ) = 34 ± 5 K and (Lir) = 7.1 x 10 12 L . There is a 
trend of increasing Td with increasing Lir for both SMGs 
and OFRGs. Firstly the SssoMm selection effectively selects 
MdXTd, thus for a given mass warmer galaxies are favoured, 
which will be more luminous. Secondly, we expect something 
of a T-L correlation because of the temperature bias in our 
luminosity estimates (see below). 

In Figure 3, we depict the locations of SMGs and 
OFRGs on the S(850^im)-z and -Lir planes. The weak cor- 
relations between S(850^tm) and Lir, (suggested by Pope 
et al. 2005) are rendered meaningless when the sample of 
OFRGs are considered together with SMGs. Selection effects 
are a concern in these diagrams, since our focus on submm- 
bright or submm-faint radio sources with relatively bright 
24^im detections, as well as SPIRE detections means that 
we may have missed populations of galaxies which would fill 
in and further elucidate these trends. 



4 DISCUSSION 

Comparing to previous studies of SMG dust temperatures, 
Chapman et al. (2005) found Td = 36 ± 7 K by using the 
FIR/radio correlation (e.g., Condon et al. 1991). As men- 
tioned earlier, these measurements relied on indirect indica- 
tors of Td . In order to more directly compare our SPIRE re- 
sults with previous studies, we applied the method of Chap- 
man et al. (2005) to our SMG sample, using only the ra- 
dio and 850/im data and fitting Dale & Helou (2002) dust 
SED templates. Td is then derived by applying a mapping 
from their Seo^m/Sioo/jm- Using the full SED information 
for these SPIRE sources, we find that Td, measured from 
5*850 and L\aghz alone, overestimates the dust tempera- 
ture by 2±2K for SMGs, and on average overestimates the 
dust temperature of OFRGs by 7±4K. This is remarkable 
agreement, given the scatter in the IR-radio relation and 
poorly defined 850/im fluxes for individual OFRGs. Kovacs 
et al. (2006) observed 21 SMGs at 350^tm and measured FIR 
luminosities independent of radio luminosity using FIR data 
ranging from observed 350/im to 1200/im. They found that 
the local FIR/radio correlation overestimates FIR luminos- 
ity by factors of ~0.2-0.4 dex for SMGs, which agrees with 
our SPIRE analysis. 

To demonstrate these discrepancies more clearly, we es- 
timate the bolometric qm, the ratio of IR integrated flux to 
radio flux, as described in Ivison et al. (2010d): 

Sir x , , SiAGHz 



qm = log( 



3.75 x 10 12 VF; 



- log( 



Wm- 2 Hi 



(2) 



where Sir taken from rest frame 8 to 1000/im as described 
before and the radio flux density Siaghz is K-corrected, 
assuming a radio slope, a — 0.75 (Ibar et al. 2010). We 
note, however, that the choice of radio spectral slope only 
affects the qm ratio by <10%. We also note that the Dale & 
Helou SED templates assume a FIR/radio correlation with 
an effective (qm)— 2.62 (Yun et al. 2001)0- Since our Sir, 
estimates are based on the greybody fits, we extrapolate 
from the peak to the 24/im point with an appropriate power- 
law, yielding (qm)= 2.43T0.2, thereby suggesting a reliable 
estimate of the infrared luminosity of z ~ 2 ULIRGs can 
be obtained from their radio luminosity using a similar q to 
that inferred locally, as inferred in Kovacs et al. (2006,2010). 

The SMG composite template (Pope et al. 2008) used 
in Fig. 1 is derived based on mid-IR to FIR data of SMGs 
to date. While it reasonably fits the SPIRE FIR data of 
the SMG sample, in some cases, it under /overestimates the 
FIR luminosities by ±0.5 dex, and unsurprisingly is much 
worse at describing the OFRG sample. This illustrates how 
a 24/im-normalized SED fitting procedure, which is common 
in the literature (Desai et al. 2009) places poor constraints 
on the breadth of FIR properties of ULIRG samples, espe- 
cially in the absence of direct FIR measurements. Similarly, 
the use of the 24/xm emission and of the Chary & Elbaz 
SED library yields an inaccurate estimate of the infrared lu- 
minosity characterized by a large scatter (~ 0.5 dex) and a 
systematic overestimation (~ x2 times) in the luminosities 



4 The Helou et al. (1988) definition of 'q' using L[40-120/mi] dif- 
fers from our adopted Ivision et al. (2010d) definition of qiR with 
L[8-1000^tm] , all other terms being the same. Thus with an aver- 
age colour conversion of 1.91, we find qiR ~ q +0.28. 
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of ULIRGs. For radio-ULIRGs under investigation here, lu- 
minosity extrapolations based on the radio emission are con- 
siderably more reliable than those based on the mid-infrared 
emission (e.g., Elbaz et al. 2010, although they have shown 
that 24^im is a good representation of Lir for lower redshift 
galaxies) . 

It is of interest to estimate the expected infrared lu- 
minosities of our SPIRE-undetected SMGs and OFRGs, al- 
though in the latter case there is some likelihood that an 
AGN is boosting the radio luminosity. The SPIRE nondetec- 
tions of SMGs are generally found to be consistent with the 
properties inferred from the detections. We use ((/ir)=2.43 
and apply the Dale & Helou SED templates normalized to 
these infrared luminosities. For the SMGs we find that for 
5 out of 11 undetected SMGs, SPIRE flux densities inferred 
using these fits lie below the detection threshold. Of the 6 
sources with predictions for SPIRE fluxes above the detec- 
tion threshold, 3 have evidence in the UV spectra for AGN 
contributions. For the OFRGs the lack of 850/im detection 
and SPIRE detections makes it likely that an AGN is gen- 
erating the radio luminosity. 

The main concern with this study is just how repre- 
sentative the sample of ULIRGS is to the wider SMG and 
OFRG populations given the use of a 24/im prior. The 24/im 
detection could favour those galaxies that are either warmer, 
or have larger fractions of warm (T>35K) gas than average, 
biasing both the luminosities and temperatures. Apart from 
misrepresenting the class of objects as a whole, the different 
temperature distributions also question the usability of lo- 
cal templates (which presumably have a more average tem- 
perature distribution). Also, because the 24/im band falls 
into a range of the rest-frame spectrum where there are a 
lot of aromatic features, it is even more complicated to as- 
sess what subclasses are selected at what redshift. While 
our study has some clear advantages compared to prior FIR 
studies of SMGs and OFRGs, the additional bias is a dis- 
advantage, motivating a more thorough investigation of the 
SPIRE properties of more complete samples. 

Our observations firstly confirm that a substantial frac- 
tion of OFRGs (-50% of the 24/mi XID SPIRE sample) 
exhibit large infrared luminosities corresponding to SFRs of 
380 M /yr (SFR [M /yr] = lxl0~ 10 L IR [L Q ], assuming 
a Chabrier IMF). The complicated selection of our sample 
in this study, including the required 24/im XID prior for 
SPIRE flux extraction, means that a precise volume den- 
sity for these IR-luminous OFRGs is not trivial to esti- 
mate. As a lower limit, given the OFRGs without 24/im 
XID prior which have yet to be studied with Herschel, the 
IR-luminous OFRGs are about 1/4 as numerous as equiva- 
lent luminosity SMGs at z ~ 2, or > 7 x 10~ 6 Mpc~ 3 for 
Lir = 2 - 6 x 10 12 L©. The SPIRE observations secondly 
confirm the higher dust temperatures for these OFRGs than 
similarly selected SMGs with SPIRE observations, consis- 
tent with our early predictions of the submm-faint radio 
populations. Understanding the different dust temperatures 
and the connections to SMGs is possible if we take the 
luminosities of SMGs and OFRGs as dominated by star- 
formation, even if an AGN contribution cannot be ruled out. 
Since L cx Md T°' J , warmer dust implies less dust. As such, 
given the similar SMG and OFRG luminosities but their dif- 
ferent typical temperatures, a possible conclusion is simply 
that OFRGs have about a third of the dust content relative 



to stellar mass. This would naturally suggest that OFRGs 
are simply a different phase of the starburst than the cooler 
SMGs. Either OFRGs are an earlier phase, before the full- 
scale build-up of dust mass, or a later phase, where much 
of the dust is blown out by superwinds or consumed by a 
central black hole. 

Our observations also clearly confirm the remarkably 
large infrared luminosities of most SMGs (Ssso/jm > 5 mjy) 
which correspond to SFRs of 710 M©/yr. The remaining 
12 SMGs in our Lockman-East catalog which are not found 
in the XID SPIRE catalog do not distinguish themselves 
in terms of 850/im or 24/im fluxes. Although their radio- 
850/im inferred Td is cooler on average than those presented 
here, it is likely source confusion in the SPIRE maps which 
limits their study. More detailed evolutionary understanding 
of high- z ULIRGs will be facilitated by future studies of deep 
Herschel observations and including the shorter wavelength 
bands. 



ACKNOWLEDGMENTS 

We thank the referee, Attila Kovacs, for a careful reading 
and constructive criticism of the manuscript. The data pre- 
sented in this paper will be released through the Herschel 
Database in Marseille HeDaM (hedam.oamp.fr/HerMES). 
SPIRE has been developed by a consortium of institutes 
led by Cardiff Univ. (UK) and including Univ. Lethbridge 
(Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. 
Padua (Italy); IAC (Spain); Stockholm Observatory (Swe- 
den); Imperial College London, RAL, UCL-MSSL, UKATC, 
Univ. Sussex (UK); Caltech, JPL, NHSC, Univ. Colorado 
(USA). This development has been supported by national 
funding agencies: CSA (Canada); NAOC (China); CEA, 
CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB 
(Sweden); STFC (UK); and NASA (USA). Data presented 
herein were obtained using the W. M. Keck Observatory, 
which is operated as a scientific partnership among Caltech, 
the University of California and NASA. The Observatory 
was made possible by the generous financial support of the 
W. M. Keck Foundation. 



REFERENCES 

Amblard, A., Cooray, A., Serra, P., Temi, P., Barton, 
E., Ncgrello, M., Auld, R., Baes, M., et al., 2010, A&A, 
aph/1005.2412 

Baugh, C. M., Lacey, C. G., Frenk, C. S., et al. 2005, MN- 
RAS, 356, 1191 

Blain, A.W., Chapman, S.C, Smail, I., Ivison R., 2004, 
ApJ, 611, 52 

Bussmann R. S., et al., 2009, ApJ, 705, 184 

Casey, C. M., Chapman, S. C, Beswick, R. J., et al. 2009, 

MNRAS, 399, 121 

Casey, C. M., Chapman, S. C, Neri, R., et al. 2010a, MN- 
RAS, arXiv0910.5756 

Casey, C. M., Chapman, S., Smail, I., et al., 2010b, MN- 
RAS, in press 

Chanial, P. et al. 2010, MNRAS, in prep 

Chapman, S. C, Richards, E. A., Lewis, G. F., Wilson, G., 

Barger, A. J., 2001, ApJ, 548, L147 



Chapman, S. C, Blain, A., Ivison, R., Smail, I., 2003a, 
Nature, 422, 695 

Chapman, S. C, Helou, C, Lewis, G. F., Dale, D., 2003b, 
ApJ, 588, 186 

Chapman, S. C, Smail, I., Blain, A. W., Ivison, R. J. 2004, 
ApJ, 614, 671 

Chapman, S. C, Blain, A. W., Smail, I., Ivison, R. 2005, 
ApJ, 622, 772 

Chapin, E., Hughes, D., Aretxaga, I., 2009, MNRAS, 393, 
653 

Chapin, E., et al, 2010 arXivl003.2647C 

Chapman, S. C, Blain, A., Smail, I., Ivison, R. J., 2005, 

ApJ, 622, 344 

Chapman, S. C, Smail, I., Blain, A., Ivison, R. J., 2004, 
ApJ, 614, 67 

Condon, J. J., Anderson, M. L., Helou, C, 1991, ApJ, 376, 
95 

Coppin, K., et al., 2006, MNRAS, 372, 1621 

Coppin, K., et al., 2008, MNRAS, 378, 214 

Daddi, E., Rottgering, H. J. A., Labbe, I., Rudnick, C, 

Franx, M., Moorwood, A. F. M., Rix, H. W., van der Werf, 

P. P., van Dokkum, P. G., 2008, ApJ, 628, 50 

Dale, D., Helou, G., Contursi, A.; Silbermann, N. A.. Kol- 

hatkar, S., 2001, ApJ, 562, 142 

Dale, D. A., Helou, G. 2002, ApJ, 576, 159 

Dave, R., Finlator, K., Oppenheimer, B. D., et al. 2010, in 

press 

Desai V., et al, 2009, ApJ, 677, 943 

Dey A., et al., 2008, ApJ, 677, 943 

Dunlop, J., et al., 2010, arXiv0910.3642D 

Eales, S., Lilly, S., Webb, T., Dunne, L., Gear, W., 

Clements, D., Yun, M., 2000, AJ, 120, 2244 

Elbaz, D., et al., 2010, A&A, in press, l|arXiv:1005.2859|l 

Greve, T., et al., 2005, ApJ, 621, 124 

Griffin, M., Abergel, A., Abreu, A., et al. 2010, A&A, in 
press 

Helou, G., Khan, I. R., Malek, L., Boehmer, L. 1988, ApJS, 
68, 151 

Ibar, E., Ivison, R., Best, P., et al., 2010 MNRAS, 401L, 
53 

Ibar, E., Ivison, R., Biggs, A., Lai, D., Best, P., Green, D. 
2009, MNRAS, 397, 281 

Ivison, R. J., et al., 2005, MNRAS, 364, 1025 

Ivison, R. J., et al., 2010a, MNRAS, 404, 198 

Ivison, R. J., et al., 2010b, MNRAS. larXiv: 1009.07491 

Ivison, R. J., et al., 2010c, MNRAS, 2010, 402, 245 

Ivison, R. J., et al., 2010d, MNRAS, 2010, arXivl005.1072 

Kovacs, A., et al. 2006, ApJ, 634, 31 

Kovacs, A., et al., 2010, ApJ, 717, 29 

Lonsdale et al. 2004, ApJS, 154, 54 

Magdis et al. 2010, MNRAS, submitted 

Nguyen et al. 2010, MNRAS, in press 

Oliver, Seb; Frost, M.; Farrah, D.; Gonzalez-Solares, E.; 
Shupe, D. L.; Henriques, B.; Roseboom, I.; Alfonso-Luis, 
A. et al. 2010, MNRAS, in press 
Pilbratt, G. et al. 2010, A&A, in press 
Pope, A., et al. 2008, ApJ, 603, L13 
Roseboom, I., et al. 2010, MNRAS, in press 
Seymour, N. et al. 2009, MNRAS, 398, 1573 
Smail, I., Ivison, R., Blain A., 1997, ApJ, 490, L5 
Wardlow, J. et al. 2010, MNRAS, submitted 



Herschel Galaxies at z>l 7 

Yun, M. S., Reddy, N. A., & Condon,J.J., 2001, ApJ, 554, 
803 

Younger J. D., et al., 2009, MNRAS, 394, 1685 



